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Abstract. We discuss three conceivable scenarios of extension and/or modification of the lAU 
relativistic resolutions on time scales and spatial coordinates beyond the Standard lAU Frame- 
work. These scenarios include: (1) the formalism of the monopole and dipole moment transforma- 
tions of the metric tensor replacing the scale transformations of time and space coordinates; (2) 
implementing the parameterized post-Newtonian formalism with two PPN parameters - /9 and 
^\ • 7; (3) embedding the post-Newtonian barycentric reference system to the Friedman-Robertson- 
' Walker cosmological model. 

o : 

_ Keywords, relativity, gravitation, standards, reference systems 

Ph! 

^ I 1. Introduction 

' Fundamental scientific program of detection of gravitational waves by the space in- 
' . terferometric gravitational wave detectors like LISA is a driving motivation for further 
I— i' systematic developing of relativistic theory of reference frames in the solar system and 
beyond. LISA will detect gravitational wave sources from all directions in the sky. These 
sources will include thousands of compact binary systems containing neutron stars, black 
holes, and white dwarfs in our own Galaxy, and merging super-massive black holes in 
distant galaxies. However, the detection and proper interpretation of the gravitational 
waves can be achieved only under the condition that all coordinate-dependent effects 
are completely understood and subtracted from the signal. This is especially important 
. for observation of gravitational waves from very distant sources located at cosmological 

■ distances because the Hubble expansion of our universe affect propagation of the waves. 
^-H ■ New generation of microarcsecond astrometry satellites: SIM and a cornerstone mission 

' of ESA - Gaia, requires a novel approach for an unambiguous interpretation of astromet- 

QQ . ric data obtained from the on-board optical instruments. SIM and Gaia complement one 

■ another. Both SIM and Gaia will approach the accuracy of 1 /las. Gaia will observe all 
0^ ■ stars 10^) between magnitude 6 and 20. The accuracy of Gaia is about 5 ^as for 

the optimal stars (magnitude between 6 and 13). SIM is going to observe 10000 stars 
with magnitude up to 20. The accuracy of SIM is expected to be a few /ias and can be 
reached for any object brighter than about 20 provided that sufficient observing time 
is allocated for that object. At this level the problem of propagation of light rays must 
' be treated with taking into account relativistic effects generated by non-static part of 
the gravitational field of the solar system and binary stars (Kopeikin & Gwinn 2000). 
Astrometric resolution in 1 yuas forces us to change the classic treatment of parallax, 
aberration, and proper motion of stars by switching to a more precise definition of ref- 
erence frames on a curved space-time manifold (Kopeikin & Schafer 1999, Kopeikin & 
Mashhoon 2000, Kopeikin et al. 2006). Advanced practical realization of an inertial ref- 
erence frame is required for unambiguous physical interpretation of the gravitomagnetic 
precession of orbits of LAGEOS satellites (Ciufolini & Pavlis 2004, Ries 2009) and GP-B 
gyroscope, which is measured relative to a binary radio star IM Pegasi that has large 
annual parallax and proper motion (Ransom et al. 2005). 
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Recent breakthroughs in technology of drag-free sateUites, clocks, lasers, optical and ra- 
dio interferometers and new demands of experimental gravitational physics (Lammerzahl 
et al. 2001, Dittus et al. 2008) make it necessary to incorporate the parameterized post- 
Newtonian formalism (Will 1993) to the procedure of construction of relativistic lo- 
cal frames around Earth and other bodies of the solar system (Klioner & Soffel 2000, 
Kopeikin & Vlasov 2004). The domain of applicability of the lAU relativistic theory of 
reference frames (Soffel et al. 2003) should be also extended outside the boundaries of 
the solar system (Kopeikin & Gwinn 2000). 

In what follows, Latin indices takes values 1,2,3, and the Greek ones run from to 
3. Repeated indices indicate the Einstein summation rule. The unit matrix is denoted 
5ij — diag(l, 1, 1) and the fully anti-symmetric symbol Cyfe is subject to £123 = 1- The 
Minkowski metric is rjap = diag(— 1, 1,1,1). Greek indices are raised and lowered with 
the Minkowski metric, Latin indices arc raised and lowered with the unit matrix. Bold 
italic letters denote spatial vectors. Dot and cross between two spatial vectors denote 
the Euclidean scalar and vector products respectively. Partial derivative with respect to 
spatial coordinates are denoted as d/dx^ or V. 



2. Standard lAU Framework 

New relativistic resolutions on reference frames and time scales in the solar system were 
adopted by the 24-th General Assembly of the lAU in 2000 (Soffel et al. 2003, Soffel 2009). 
The resolutions are based on the first post-Newtonian approximation of general relativity. 
They abandoned the Newtonian paradigm of space and time and required corresponding 
change in the conceptual basis and terminology of the fundamental astronomy (Capitainc 
et al. 2006). 

Barycentric Celestial Reference System (BCRS), x" = (ct,a;), is defined in terms of a 
metric tensor gap with components 

5oo = -l + ^-^ + 0(c-^), (2.1) 
5o, = -^ + 0(c-5), (2.2) 

5..=%(l + 5^)+0(c-4). (2.3) 

Here, the post-Newtonian gravitational potential w generalizes the Newtonian potential, 
and is a vector potential related to the gravitomagnetic effects (Ciufolini & Wheeler 
1995). These potentials are defined by solving the field equations 

Uw = -AnGa , (2.4) 
Dw' = -47rGcr* , (2.5) 

where □ = -c-'^d'^/dt^ + is the wave operator, a = c-'^{T°° + T'*"), cr* = c'^T^', 
and T'"' are the components of the stress-energy tensor of the solar system bodies. 

Equations (2.4), (2.5) are solved by iterations 

^{t,x)=Gj j Sx'a{t,x')\x-x'\+0{c-^), (2.6) 

w\t, x)^G f + 0{cr') , (2.7) 
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which are to be substituted to the metric tensor (2.1)-(2.3). Each of the potentials, w 
and w', can be Hnearly decomposed in two parts 

w = ive + w , (2-8) 
w' ^w'e+w\ (2.9) 

where We and We BCRS potentials depending on the distribution of mass and current 
only inside the Earth, and we and We gravitational potentials of external bodies. 

Geocentric Celestial Reference System (GCRS) is denoted X" = (cT,X). OCRS is 
defined in terms of the metric tensor Gap with components 

Goo = -l + -^-^ + 0{c-') , (2.10) 

Go. = --^ + 0(c-5), (2.11) 

G,, = % (^1 + ^) + 0(c-4) . (2.12) 

Here W = W{T,X) is the post-Newtonian gravitational potential and W^{T,X) is a 
vector-potential both expressed in the geocentric coordinates. They satisfy to the same 
type of the wave equations (2.4), (2.5). 

The geocentric potentials are split into two parts: potentials We and We arising from 
the gravitational field of the Earth and external parts associated with tidal and inertial 
effects. lAU resolutions implied that the external parts must vanish at the geocenter and 
admit an expansion in powers of X (Soffel et al. 2003, Soffel 2009) 

w^(r,x) = VFB(T,x) + Wki„(T,x) + w^dyn(r,x) , (2.13) 

W\T,X) = Wi,{T,X) + VF^i„(T,X) + W^y^{T,X) . (2.14) 

Gcopotentials We and arc defined in the same way as we and w'e (see equations 
(2.6)--(2.7)) but with quantities cr and cr' calculated in the GCRS. Wkin and W^^^^ are 
kinematic contributions that arc linear in spatial coordinates X 

1 
4 

where Qi characterizes the minute deviation of the actual world line of the geocenter 
from that of a fiducial test particle being in a free fall in the external gravitational field 
of the solar system bodies (Kopeikin 1988) 

Q,=dMxE)'a'E + Oic-^) . (2.16) 

Here = dv\^/dt is the barycentric acceleration of the geocenter. Function Op^^^. de- 
scribes the relativistic precession of dynamically non-rotating spatial axes of GCRS with 
respect to reference quasars 

The three terms on the right-hand side of this equation represent the geodetic, Lcnse- 
Thirring, and Thomas precessions, respectively (Kopeikin 1988, Soffel et al. 2003). 

Potentials Wdyn and VFjyjj are generalizations of the Newtonian tidal potential. For 
example, 

Wdy„(r, X) = w{xE + X) - w{xe) - X'dMxE) + 0{c-^). (2.18) 
It is easy to check out that a Taylor expansion of w{xe + X) around the point xe yields 



w^,^ = Q,x\ wi,^ = -c\,p,{np -ni^,,)x\ (2.15) 
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VVdyn(T, X) in the form of a polynomial starting from the quadratic with respect to X 
terms. We also note that the local gravitational potentials We and of the Earth are 
related to the barycentric gravitational potentials we and by the post-Newtonian 
transformations (Brumberg & Kopeikin 1989, Soffcl et al. 2003). 



3. lAU Scaling Rules and the Metric Tensor 

The coordinate transformations between the BCRS and GCRS are found by matching 
the BCRS and GCRS metric tensors in the vicinity of the world line of the Earth by 
making use of their tensor properties. The transformations are written as (Kopeikin 
1988, Soffel et al. 2003) 



T = 



1 



ve ■ rE] + ^ 



B'^r'Er'E + C'{t,x) 



-VeVe ■ ve + w{xE)rE + '^eO'E ■ rE 



'■E' E 



+ 0(c-5), 
0(c-4), 



where te = x — xe, and functions A, B, , C{t, x) are 



(3.1) 
(3.2) 

(3.3) 
(3.4) 
(3.5) 
(3.6) 
(3.7) 



Here again x^, v^^, and arc the barycentric position, velocity and acceleration vectors 
of the Earth, the dot stands for the total derivative with respect to t. The harmonic 
gauge condition does not fix the function C{t, x) uniquely. However, it is reasonable to 
fix it in the time transformation for practical reasons (Soffel et al. 2003). 

Earth's orbit in BCRS is almost circular. This makes the right side of equation (3.3) 
almost constant with small periodic oscillations 



dA 1 


v'e + w{xe), 






'dt ~ 2 






dB 

m ^ ^ 


1 4 3 -/ \ 

g^E~ 2^eW[xe) 


+ 4v'e w' + ^ 


w^{xe), 


B' = - 


^vlv^E+iw'ixE) 


- 3vEwixE) 


+ ]^5''^w{xe), 


B'^ = - 


v'eQj +2djw\xE) 


- VEdjW{xE) 


C{t,x) = - 


— rE (aE - rE). 







— v%'{-w{xe)~ c?Lc + (periodic terms) 



(3.8) 



where the constant Lc and the periodic terms have been calculated with a great precision 
in (Fukushima 1995). For practical reason of calculation of ephemerides of the solar 
system bodies, the BCRS time coordinate was re-scaled to remove the constant Lq from 
the right side of equation (3.8). The new time scale was called TDB 



troB — i (1 ~ Lb) , 



(3.9) 



where a new constant is used for practical purposes instead of Lc in order to take into 
account the additional linear drift between the geocentric time T and the atomic time 
on geoid, as explained in (Brumberg & Kopeikin 1990, Irwin & Fukushima 1999, Klioncr 
et al. 2009). Time re-scaling changes the Newtonian equations of motion. In order to 
keep the equations of motion invariant scientists doing the ephemerides also re-scaled 
spatial coordinates and masses of the solar system bodies. These scaling transformations 
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are included to lAU 2000 resolutions (SofFel et al. 2003) but they have never been explic- 
itly associated with transformation of the metric tensor. It had led to a long-standing 
discussion about the units of measurement of time, space, and mass in astronomical 
measurements (Klioner et al. 2009). 

The scaling of time and space coordinates is associated with a particular choice of the 
metric tensor corresponding to either TCB or TDB time. In order to see it, we notice 
that equation (2.15) is a solution of the Laplace equation, which is defined up to an 
arbitrary function of time Q. If one takes it into account, equations (2.15), (3.3) can be 
re-written as 

Wkin = Q + Q^X' , (3.10) 

^-lvl+wixE)-Q, (3.11) 

and, if we chose Q — c^Lc, it eliminates the secular drift between times T and t with- 
out explicit re-scaling of time, which is always measured in SI units. It turns out that 
Blanchet-Damour (Blanchet &; Damour 1989) relativistic definition of mass depends on 
function Q and is re-scaled in such a way that the Newtonian equations of motion remain 
invariant. Introduction of Q to function W brings about implicitly re-scales spatial co- 
ordinates as well. We conclude that introducing function Q = <?Lc to the metric tensor 
without apparent re-scaling of coordinates and masses might be more preferable for lAU 
resolutions as it allows us to keep the SI system of units without changing coordinates 
and masses made ad hoc " by hands" . Similar procedure can be developed for re-scaling 
the geocentric time T to take into account the linear drift existing between this time and 
the atomic clocks on geoid (Brumberg & Kopeikin 1990, Irwin & Fukushima 1999). 



4. Parameterized Coordinate Transformations 

This section discusses how to incorporate the parameterized post-Newtonian (PPN) 
formalism (Will 1993) to the lAU resolutions. This extends applicability of the resolu- 
tions [ to a more general class of gravity theories. Furthermore, it makes the lAU resolu- 
tions fully compatible with JPL equations of motion used for calculation of ephemerides 
of major planets. Sun and Moon. 

These equations of motion depend on two PPN parameters, (3 and 7 (Seidelmann 
1992) and they are presently compatible with the lAU resolutions only in the case of 
/? = 7 = 1. Rapidly growing precision of optical and radio astronomical observations as 
well as calculation of relativistic equations of motion in gravitational wave astronomy 
urgently demands to work out a PPN theory of relativistic transformations between the 
local and global coordinate systems. 

PPN parameters (3 and 7 are characteristics of a scalar field which makes the metric 
tensor different from general relativity. In order to extend the lAU 2000 theory of ref- 
erence frames to the PPN formalism we employ a general class of Brans-Dicke theories 
(Brans & Dicke 1961). This class is based on the metric tensor g^fj and a scalar field (f> 
that couples with the metric tensor through function 6{(j)). We assume that (f) and 6'(0) 
are analytic functions which can be expanded in a Taylor series about their background 
values (f> and 6. 

The parameterized theory of relativistic reference frames in the solar system is built 
in accordance to the same rules as used in the lAU resolutions. The entire procedure is 
described in our paper (Kopeikin & Vlasov 2004). The parameterized coordinate trans- 
formations between BCRS and GCRS arc found by matching the BCRS and GCRS 
metric tensors and the scalar field in the vicinity of the world line of the Earth. The 
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transformations have the foUowing form (Kopeikin & Vlasov 2004) 

T = t-\[A + VE-rE] + ^\B + BW'E + B'^ r'^ t'e + C{t, x)] + Oicr^), (4.1) 



' E 



^'We^E + iQ^E + IwixEVE + rElE^E 



E' E 



where ve = x — xe, and functions A{t),B{t),B^{t),B'^^ {t),C{t,x) are 
dA 



^v%+w-Q{xe), 



+ 0(c-4)(4.2) 



(4.3) 



dB 
'dt 



1 



1 



-5r = -o4- 7+7T vj,wixE) + 2{l+j)v},w'+{P-- w'ixE), (4.4) 



1 



1 



-§ 4 «B + 2(1 + 7) w'ixE) - (1 + 27) v}, w{xe), 



B' 

B'^ = -VeQj + {l+'^)djW'{xE) -jv'Ed.jw{xE) + ^S'^w{xe), 



1 



C{t,x) 



-^4 (as • ^e) 



(4.5) 
(4.6) 
(4.7) 



These transformations depends expUcitly on the PPN parameters (3 and 7 and the scahng 
function Q, and should be compared with those (3.1)-(3.7) adopted in the lAU resolu- 
tions. 

PPN parameters /3 and 7 have a fundamental physical meaning in the scalar-tensor 
theory of gravity along with the universal gravitational constant G and the fundamental 
speed c. It means that if the parameterized transformations (4.1)-(4.7) are adopted by the 
lAU, the parameters /? and 7 must be included to the number of the astronomical con- 
stants which values must be determined experimentally. The program of the experimental 
determination of P and 7 began long time ago and it makes use of various observational 
techniques. So far, the experimental values of /? and 7 are indistinguishable from their 
general-relativistic values (3 — 1,^—1. 



5. Matching lAU Resolutions with Cosmology 

BCRS assumes that the solar system is isolated and space-time is asymptotically flat. 
This idealization will not work at some level of accuracy of astronomical observations 
because the universe is expanding and its space-time is described by the Friedman- 
Robertson- Walker (FRW) metric tensor having non-zero Ricmannian curvature (Misner 
et al. 1973). It may turn out that some, yet unexplained anomalies in the orbital motion of 
the solar system bodies (Anderson 2009, Krasinsky & Brumberg 2004, Lammerzahl et al. 
2006) are indeed associated with the cosmological expansion. Moreover, astronomical 
observations of cosmic microwave background radiation and other cosmological effects 
requires clear understanding of how the solar system is embedded to the cosmological 
model. Therefore, it is reasonable to incorporate the cosmological expansion of space- 
time to the lAU 2000 theory of reference frames in the solar system (Kopeikin et aL2001, 
Ramirez et al. 2002). 

Because the universe is not asymptotically-flat the gravitational field of the solar sys- 
tem can not vanish at infinity. Instead, it must match with the cosmological metric tensor. 
This imposes the cosmological boundary condition. The cosmological model is not unique 
has a number of free parameters depending on the amount of visible and dark matter, 
and on the presence of dark energy. We considered a FRW universe driven by a scalar 
field imitating the dark energy and having a spatial curvature equal to zero (Dolgov 
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et al. 1990, Mukhanov 2005). The universe is perturbed by a localized distribution of 
matter of the solar system. In this model the perturbed metric tensor reads 

gap = a^{v)faP , /a/3 = {Vap + hap) , (5.1) 

where perturbation hap of the background FRW metric tensor g^p ~ a^"(]ap is caused by 
the presence of the solar system, a(r]) is a scale factor of the universe depending on the, 
so-called, conformal time 77 related to coordinate time t by simple differential equation 
dt = a{r])dr]. In what follows, a linear combination of the metric perturbations 

^ap ^ f^ap _ Iri^l^h , (5.2) 

where h = rj^^hap, is more convenient for calculations. 

We discovered a new cosmological gauge, which has a number of remarkable proper- 
ties. In case of the background FRW universe with dust equation of state (that is, the 
background pressure of matter is zero (Dolgov et al. 1990, Mukhanov 2005) this gauge 
is given by (Kopcikin et a/.2001, Ramirez et al. 2002) 

r^ip - 2H^5^ , (5.3) 

where bar denotes a covariant derivative with respect to the background metric gap, 
(fi = (/)/a^, H = d/a is the Hubble parameter, and the overdot denotes a time derivative 
with respect to time r]. The gauge (5.3) generalizes the harmonic gauge of asymptotically- 
flat space-time for the case of the expanding non-flat background universe. 

The gauge (5.3) drastically simplifies the linearized Einstein equations. Introducing 
notations 700 = 4w/c^, 70; = —iw^/c^, and 7^ = Aw^^ /c'^, and sphtting Einstein's 
equations in components, yield 



□X - 2Hd^x + \h^X = -4^Ga , (5.4) 

Uw - 2Hdr,w = -47rGtT - AH^x , (5.5) 

□w* - 2Hdrjw' + H^w' = -AirGa' , (5.6) 

Ow'^ - 2Hd„w'^ = -AttGT^ , (5.7) 



where 9^ = d/drj, □ = — c + V^, x = w — <p/2, the Hubble parameter H = d/a = 
2/r], densities <t = c-2(T"" + T"'), a' ^ c'^T^' with T°'^ being the tensor of energy- 
momentum of matter of the solar system defined with respect to the metric fap. These 
equations extend the domain of applicability of equations (2.4), (2.5) of the lAU standard 
framework (Soffcl 2009) to the case of expanding universe. 

First equation (5.4) describes evolution of the scalar field while the second equation 
(5.5) describes evolution of the scalar perturbation of the metric tensor. Equation (5.6) 
yields evolution of vector perturbations of the metric tensor, and equation (5.7) describes 
generation and propagation of gravitational waves by the isolated N-body system. Equa- 
tions (5.4)-(5.7) contain all corrections depending on the Hubble parameter and can be 
solved analytically in terms of generalized retarded solution. Exact Green functions for 
these equations have been found in (Kopcikin et aL2001, Ramirez et al. 2002, Haas & 
Poisson 2005). They revealed that the gravitational perturbations of the isolated system 
on expanding background depend not only on the value of the source taken on the past 
null cone but also on the value of the gravitational field inside the past null cone. 

Existence of extra terms in the solutions of equations (5.4)-(5.7) depending on the 
Hubble parameter brings about cosmological corrections to the Newtonian law of grav- 
ity. For example, the post-Newtonian solutions of equations (5.5), (5.6) with a linear 
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correction due to the Hubble expansion are 

Hv,-) = Gj + (5-8) 

-GH J d^x'cr{7], x') + O {c-'^) + O (H^) , 

(77, x) = G [ '"'i'^'^'^f'^' + O (c-2) + O (H^) . (5.9) 
J \x ^ x'\ 

Matching these solutions with those defined in the BCRS of the lAU 2000 framework 
in equations (2.6), (2.7) is achieved after expanding all quantities depending on the 
conformal time rj in the neighborhood of the present epoch in powers of the Hubble 
parameter. 

Current lAU 2000 paradigm assumes that the asymptotically-flat metric, fap, is used 
for calculation of light propagation and ephemerides of the solar system bodies. It means 
that the conformal time rj is implicitly interpreted as TCB in equations of motion of light 
and planets. However, the physical metric gafs differs from /^^ by the scale factor a^ir]), 
and the time rj relates to TCB as a Taylor series that can be obtained after expanding 
a{ri) = flo + 077 + in polynomial around the initial epoch ijq and defining TCB at the 
epoch as TCB=ao77. Integrating equation dt = a{r])dr] where t is the coordinate time, 
yields 

t^TGB + ^H-TCB^ + ... , (5.10) 

where H = H/qq and ellipses denote terms of higher order in the Hubble constant Ti. 
The coordinate time t relates to the atomic time TAI (the proper time of observer) 
by equation, which does not involve the scale factor 0(77) at the main approximation. 
It means that in order to incorporate the cosmological expansion to the equations of 
motion, one must replace TCB to the quadratic form 

TCB — > TCB + • TCB^ . (5.11) 

Distances in the solar system are measured by radio ranging spacecrafts and planets. 
Equations of light propagation preserve their form if one keeps the speed of light con- 
stant and replace coordinates {ri,x) to {t,a), where the spatial coordinates H relate to 
coordinates x by equation da — a{i])dx. Because one uses TAI for measuring time, the 
values of the spatial coordinates in the range measurements are given in terms of the 
capitalized coordinates H. Therefore, the ranging measurements arc not affected by the 
time transformation (5.10) in contrast to the measurement of the Doppler shift, which 
deals with time only. This interpretation may reconcile the " Quadratic Time Augmenta- 
tion" model of the Pioneer anomaly discussed by Anderson et al (2002) in equation (61) 
of their paper. 
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